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Abstract. Let (p n )n be a sequence of orthogonal polynomials with respect 
to the measure fi. Let T be a linear operator acting in the linear space of poly- 
nomials P and satisfying that deg(T(p)) = deg(p) — 1, for all polynomial p. We 
then construct a sequence of polynomials (s n ) n , depending on T but not on fi, 
such that the Wronskian type n X n determinant det (T I— 1 (p m +j_i (a;))) ? 

(—1 \ m 

qL , up to multiplica- 
n+I / i,j=l 

tive constants, where the polynomials q l n , n,i > 0, are defined by q l n (x) = 
~^2"—q lA Sn—j (s), and fj,j are certain generalized moments of the measure fi. 
For T = d/dx we recover a known result by Leclerc which extends the well- 
known Karlin and Szego identities for Hankel determinants whose entries are 
ultraspherical, Laguerre and Hermite polynomials. We also show that for cer- 
tain operators T, the second determinant above can be rewritten in terms of 
Selberg type integrals, and that for certain operators T and certain families 
of orthogonal polynomials (p n ) n , one (or both) of these determinants can also 
be rewritten as the constant term of certain multivariate Laurent expansions. 
As a consequence, we get some new Selberg type integrals and constant terms 
identities (as well as new proofs for classical ones). 



1. Introduction and results 

Determinants whose entries are orthogonal polynomials is a long studied subject. 
One can mention Turan inequality for Legendre polynomials |21] and its generaliza- 
tions, specially that of Karlin and Szego on Hankel determinants whose entries are 
ultraspherical, Laguerre, Hermite, Charlier, Meixner, Krawtchouk and other fam- 
ilies of orthogonal polynomials ([II])- Karlin and Szego's strategy was to express 
these Hankel determinants in terms of the Wronskian of certain orthogonal poly- 
nomials of another class. However the proofs of Karlin and Szego did not clearly 
show what in their formulas resulted from a general algebraic transformation, and 
what in contrast was due to some particular properties of the orthogonal polyno- 
mials under consideration. Leclerc ([H]) clarified this by stating that a Wronskian 
of orthogonal polynomials is proportional to a Hankel determinant whose elements 
form a new sequence of polynomials (not necessarily orthogonal). More precisely, 
let (a n ) n be a sequence of numbers. With this sequence Leclerc associated two 
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classes of polynomials, p n , q n , n > 0, denned by 



(1.1) 



(1.2) 



«0 
«1 



O-n-1 
«.,, 



/ j 3 
3=0 



a i 

a 2 

n 



ft,, 



&1n-1 



1 

X 

„n-l 



(- x y 



It is well- known that if a n , n > 0, are the moments of a non-degenerate positive 
measure fi, a n = Jt n d[i, the polynomials (p n ) n are then orthogonal with respect 
to fJi. 

Theorem 1.1. (Theorem 1 of [15] ) The following identity holds for all integer 
values of n and m, n,m > 1, 

(1-3) dct = C„, m det (q n+i+j _ 2 (x))™ j=1 , 

where C n , m is independent of x: 



{-l) nm \{k\Aet{a l+ ^ 2 ) k i 



k-\-m 



fc=i 



J.j=n X 3 



1 if m < n and det(C) = 1 if the size of C reduces 



As usual, we take ]J j 
to x 0. 

In 0], the author has conjectured some identities for Casorati determinants of 
classical discrete orthogonal polynomials. These conjectures can be considered as 
a discrete version of the Karlin and Szego's identities for Hankel and Wronskian 
determinants of ultraspherical, Laguerre and Hcrmitc polynomials. 

Here it is one of these conjectures. Let (c^)„, a ^ 0, be the monic Charlier 
polynomials (see (|5.14p below), then for n, k, m > and a £ R \ {0}, we have 



(1.4) 



det (Sg+i-ifa + i - . =1 _ det {t-^i-k + i - 1)) 



In fact, Karlin and Szego considered in their paper |11) the determinant in the left 
hand side of the previous formula (and the corresponding for Meixner, Krawtchouk 
and Chebyshev) , but they were more interested in studying the sign of this deter- 
minant that in computing it explicitly. 

The purpose of this paper is to show that Leclerc identities for Wronskian of or- 
thogonal polynomials and author conjectures for Casorati determinants of classical 
discrete orthogonal polynomials are particular cases of a quite more general result. 

Indeed, consider a linear operator T acting in the linear space of polynomials P 
and satisfying that deg(T(p)) = deg(p) — 1, for all polynomial p. We associate two 
sequences of polynomials (r„)„ and {s n ) n to this operator T. 

(1) The sequence (r n ) n satisfies that 

(a) ro = 1 and the degree of r n is n, n > 0, 

(b) T(r n ) = r n _i, n > 1. 
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It is easy to check that this sequence of polynomials always exists and it is 
unique if we fix the values of r n , n > 1, at a given number xo (see Section 
4 for more details). 
(2) The sequence (s n ) n are now defined recursively by sq = 1 and 

n 

(1.5) j2 s ^ x >«-jW = °- 

3=0 

We are now ready to establish the main result of this paper. 

Theorem 1.2. Consider a linear operator T acting in the linear space of poly- 
nomials P and satisfying that deg(T(p)) = deg(p) — 1, for all polynomial p. We 
associate it the two sequences of polynomials (r n )„ and (s n ) n as above. Let [i be a 
measure and consider a sequence (p n ) n of orthogonal polynomials with respect to fi 
(see Section 2 for more details). For a given sequence of polynomials (tpiji, ipi of 
degree i, we write (£, i,j > 0, for the numbers fj,j — J rjipidpL. We now consider 
the polynomials q\, i,n > defined by 

n 

(1.6) q i n (x)=^n i j s n -j(x). 

3=0 

We then have 

(1.7) fi m _ 1 det(T i - 1 (p m+j _ 1 (aO)) 4 B =C„, m det (ct+i-x^Y > 
where fl m —i and C n . m are independent of x: 

n— 1 

n m _x = det( M ^)™ =1 , C n , m = (-!)-+(?) J] ^±1 

j=Q cr ro+j 

and and <j n are the coefficients of x n in the power expansion of p n and r n , 
respectively. 

Let us note that we have some degrees of freedom in the polynomials q l n , n, i > 0. 
Indeed, for a given number xq, they depend on the numbers r n (xo), n > 1, and on 
the sequence of polynomials (ipi)i. 

The content of this paper is as follows. In Section 4 we prove Theorem 11.21 
(actually a further extension for orthogonal polynomials with respect to a bilinear 
form defined in P x F). The proof is similar to that of Leclerc for Theorem 11.11 
and it is based on some Turnbull determinantal identities (see Section 3). Section 
5 is devoted to specialize Theorem 11.21 to some relevant operators (such as the 
derivative, the first order difference operator A, the q-derivative D q , and others). 
We show that for these operators, the polynomials q % n and the determinant in the 
right hand side of (|1.7[) have a richer estructure and more properties. In particular, 
the specialization for the first order difference operator A will allow us to prove our 
conjectures in [3] for Casorati determinants of Charlier, Meixner and Krawtchouk 
orthogonal polynomials. 

In Section 6 we show that for certain operators T, the determinant in the right 
hand side of the previous theorem can be rewritten in terms of Selberg type inte- 
grals. Indeed, write A(x) for the Vandermonde determinant 

(1.8) A(x)= [J (xi-xj), 

l<i<j<m 
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where x stands for the multivariable x = (xi, • • • , x m ). Given a measure /x, consider 
the multiple integral 

(1.9) /-/^(xWx), 



where cfyx(x) = cfyx(xi) • • ■ d)j,(x m ) and s(x) is a symmetric polynomial in the vari- 
ables x±, ■ ■ ■ ,x m . It seems to be rather difficult to get explicit fomulae for the 
integral (jl.9[) even for simple measures fj,. 

When /x is the Jacobi measure d\i — x a ~ l (\ — x) l3 ~ 1 dx, a, (3 > 0, and s(x) = 1, 
the integral (jl.9p is the celebrated Selberg integral ([19], see also [7]): 



Jo Jo [ 3 1 j) H r(a + /3 + (m + i-l)7)T(l + 7) ' 

Other choices of s(x) give well-known extensions of the Selberg integral. For in- 
stance, s(x) = Y\^=i{xj— u) is Aomoto's integral pQ, ands(x) = Ja(xi, • • ■ ,x m , 1/7), 
where 7 € N and J\ is a Jack polynomial, is Kadell's integral [9 . 

For s(x) = 1 and fj, the Hermite measure dpi = e~ x / 2 dx, the integral (|1.9[) is the 
also celebrated Mehta integral. 

Given a number u and a linear operator T as in Theorem 11.21 we associate it 
the two sequences of polynomials (r n ) n and (s n ) n as above. In this case we take 
T n {u) = 0, n > 1 (these conditions define r n , n > 0, uniquely). We write r„ = r n ,u 
to stress the dependence of r„ on u. We then have. 

Theorem 1.3. Assume that there are numbers aj, j > 1, with a± = u and r n u (x) — 
a n IljLi ( x ~ a j)' n>\. Let (p n ) n be the sequence of monic orthogonal polynomials 
with respect to the measure fi, then 
p p m 

(1-10) J ■ J A 2 (x) Y[r»A x j)Mx) = C T ,n, m det (T J - 1 (p m+J _ 1 ( u )))" j=1 , 



where 



c T ,n, n - (-i)- m m!<» n iip.ii 2 n ^ 

3=0 3=0 



Tn — 1 n — 1 



For n = 1, the identity (flTTO]) is known (see [20], (2.2.10)). 

Theorem 11.31 applies to the three relevant operators T = d/dx, T — A and 
T = Dq since in each case we can choose polynomials r„, n > 0, with T(r n ) = r„_i 
and satisfying the hypothesis in Theorem 11.31 We can then deduce some Selberg 
type integrals and sums corresponding to 7 = 1 in (|1.9|) . For instance, for T = d/dx, 
n = 0, 11 = 0, and /i the Jacobi weight, we recover the Selberg integral for 7=1. The 
case n = 1 and again 7 = 1 is Aomoto's result for (|1.9I) when s(x) = njlil^j ~ u ) 
(PQ). For T = A, n — 0, u = and /i the Meixner, Krawtchouk or Hahn weights 
Theorem 11.31 gives a deduction for Askey conjectures 5, 6 and 7, respectively, in [2] 
when 7 = 1. The case n = 1 and again 7=1 can be considered a discrete version 
of Aomoto's integral. Similar results can be obtained for some of the q-classical 
weights and Askey conjectures for q-Selberg integrals. 

Theorem 11.31 also gives the following Selberg type sum related to the Racah 
weight (as far as the author knows it has not been considered previously in the 
literature). Let N,j3,-y,6 be real numbers satisfying that N is a positive integer 
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and 5, 7, -/3 + 7 ^ -1, • • • , -N, ft + S ^ -1, • • • , -27V - 1. Then for m > 1 and 
< n < N 

A2/ xTT (-iV, ,8 + 5 + 1,7 + 1,(7 + ^ + 3)72),. (7 + 5 + l),. +n 



^ (JV + 7 + *+2,- i 8 + 7+l,(7 + * + l)/2 ) *+l) s ,(a: i -n)! 



= f _ ir „ r tt (j + l)!(^-7-^,-^-^,/8 + l)j(-^,/3 + ^+l,7 + l)»+j 
1 ] f} W-N+lfotf-N + fijtf-N + m + j)* 

where x A = (X(xi), • • • , \{x m )), \{x) = X ( X + 1 + S + 1), M = 

{ux, • • • , Ufe)j = (ui)j ■ ■ ■ (uk)j and (a) n stands for the Pochhammer symbol defined 

by (a)o = 1, (a) n = a(a + 1) • • • (a + n — 1), n > 1. 



In Section 7, we show that for certain operators T and certain families of orthog- 
onal polynomials (p n ) n , one (or both) of the determinants in (|1.7[) can be rewritten 
as the constant term of certain multivariate Laurent expansions. This is not sur- 
prising, because certain identities for constant term of Laurent polynomials are 
closely related to Selberg integrals. The first of these constant term identities was 
conjectured by Dyson [5] in 1962 in his study of the statistical theory of complex 
systems: 

where we denote by C.T.| z= o-F(z) the constant term of the Laurent expansion at 
z = («!,... , Zn ) = ((),■•• ,0) of F(z). 

Macdonald [16] in 1982 related Dyson identity with the root system A n -\ and 
posed some other conjectures which generalize Dyson constant term identity to 
all finite root systems. Morris |17] also proposed some constant term identities 
related to root systems and Selberg integrals. The following one can be considered 
a contour version of Selberg integral: 

(1-11) 

C - T -l-0 JUL z rn +H n-l) A W ^> 11 ( x - m + kj)\(m + kj)\k\' 

where x,m,k,n are nonnegative integers, m < x and 
For T — A and /i the Meixner measure 

^ a ^r(a; + c) . 

^ = 2^ — -} — **> 

Theorem 11.21 gives a generalization of Morris identity (|l.lll) for k = 1. Indeed, 
using the integral representation for the Meixner polynomials (m^' c )„ provided by 
the generating function 

oo 

--T (l-z)- x - c = Ym^{x)z\ 



n=0 
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we prove the following constant term identity: 

\2 



(i--r 



CT.L-oTT — ^-r TT U 

3 > j l<i<j<n 

= (_i)(?)+G)+»™ n ! a (?) det (m£?ri n-i ~ i+3 (-x + j - I))™ . 

V ) i,j=l 

For c = —x — n + 1, the determinant in the right hand side of this identity can be 
easily computed, and we recover Morris identity (jl.lip for k = 1. For m = 1, we 

get 

" (i _ ^jl) x 

c - T -i-°n (1 - .^i-^n n (^-^) 2 =(-i)^n!-r c -" +i (^)- 

j=l V Jj j l<i<j<n 

which it can be considered a contour integral version of Aomoto's integral pQ. 

We find other constant term identities using ultraspherical and Charlier polyno- 
mials. 



2. Preliminaries 



Let (, } be a bilinear form acting in the linear space of polynomials: (, ) : P x 
P — > C. Let K be a positive integer or infinity. We say that the polynomials p n , 
< n < K, are pseudo-orthogonal with respect to (, } if (p n ,pj) = 0, < j < n — 1 
and 1 < n < K . If, in addition, p n has degree n, < n < K, we say that (p n )n=o is 
a sequence of orthogonal polynomials with respect to (, ). It is clear that the n-th 
orthogonal polynomial with respect to a bilinear form, if there exists, is unique up 
to multiplicative constants. 

Let (r„)^L and (tpn)n=o be two sequences of polynomials satisfying that dcg(r„) = 
deg(V>n) = n. Write a n for the leading coefficient of r n . Write also f/j — (rj,ipi), 
< i,j < K. We can then construct pseudo orthogonal polynomials p n , n — 
0, • ■ ■ , K, with respect to (, ) using the formula 



(2.1) Pn (x) 



11° 
11° 



/'() 



Mn-1 



/'o 



r n {x) 

1 



1 < n < K, po = l. 



Indeed, from the definition, we straightforwardly have that (p n ,ipj) = 0, j = 
0, • ■ ■ , n — 1. Since deg(rj) = j, we have that deg(p n ) < n, and since deg^) = j 
as well, we can write each pj as a linear combination of tpi, i = 0, • ■ • , j. As a 
consequence we get (p n ,Pj) = 0, j = 0, • • ■ ,n — 1. 

Consider now the determinants Q n , n > 0, defined by 



(2.2) 



/4! 



It is clear from (|2.ip that the coefficient of x n in the power expansion of p n is 
(-l)"O n _ 1( r n . 

Hence, the bilinear form (, ) has a sequence p n , n = 0, ••■ ,K, of orthogonal 
polynomials if and only if f2„ ^ 0, < n < K. We can then construct a sequence 
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of monic orthogonal polynomials (p n )£L , p n of degree n, using the formula 
(2-3) 



Pn{x) 



(-1)" 
^ n — 1 &n 



M° Mn 1 r n (x) 

,0 ,,1 ..n— 1 



Mo 



/'() 



MS" 1 



1 



1 < n < if, po = 1- 



The determinant fi„ can be computed in terms of the norm of the monic poly- 
nomials. Indeed, write ||j>n|| 2 = (pmPn)- If we write v n for the leading coefficient 



of tpn, we have ip r , 



V n Pr, 



YJj=o a njPj- And hence (p n ,ip n ) 



Using 



, we also get (p n ,ip n ) = (— l) n f2 n /(Q n _i(7 n ). From where we get 



(2.4) 



(-1) 



n(ri+l)/2 



lift- 1 



In this paper, we are mainly interested in bilinear form defined by a measure \x 
supported in the real line: {p, q) = J pqd^i. 

A bilinear form can be represented by a measure supported in the real line if and 
only if it is Hermitian ((p, q) = (q,p), p, q G P) and the operator of multiplication 
by x is symmetric with respect to the bilinear form: (xp, q) — (p, xq) , p, q € P. 

We need some combinatorial identities which we include in the following lemma. 
Lemma 2.1. For x G K, n, jeN with i < g < n + i, we have 



(2.5) 
(2.6) 



—x + i 
n + i — g 



min(n,p) 

E 



J 

I - 1 



x-j 

I 



i + j - 9 



9 

9- j 



—x 
n-j 



Proof. Both identities can be easily proved from the basic identity (") = ( t ' t , 1 ) + 
("Zi) using induction on x. 



□ 



3. LECLERC VERSION OF TURNBULL DETERMINANTAL IDENTITIES 

In this section, we recall briefly Leclerc notation of 14 for minor identities (which 
is a variant of Turnbull's dot notation, [25], p. 27). Let M be a n x p matrix with 
p > n and a, 6, c, • • • , d be n column vectors of M. The maximal minor of M taken 
on these n columns is denoted by either a bracket or a one line tableau: 



[ab- 



A product of k minors of M is designated by a k x n tableau. For instance, for 
k = 2: 



a 6 


c 


e 


• / 



To denote alternating sums of products of minors, tableaux with boxes enclosing 
certain vectors are used. Let T be a k x n tableau and A a subset of elements of T. 
Given a permutation a of the elements of A, denote by a(T) the tableau in which 
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the elements of A are permuted by a. Now, boxing in T the elements of A, we 
get a new tableau r which will represent the alternating sum of all tableaux cr(T), 
taking into account the fact that a permutation of elements of the same row gives 
a trivial action. More precisely r is defined by 



r = ^sign( ( r)(7(T) I 



where a runs trough the cosets of the symmetric group S(A) modulo the subgroup 
of the permutations which leave unchanged the rows of T. For example, 






a b c 




/ b 


c 




a f c 




d e / 




d e 


a 




deb 



We can now state Turnbull's identity (see [22], p. 48 and [14]). 

Theorem 3.1. Letr be apxn tableau with set of enclosed elements A of cardinality 
< n. Let R be one of the rows of T and denote by B the set of elements of R which 
are not enclosed. Form a new tableau v by (i) exchanging the elements of A which 
do not belong to R with elements of B; (ii) removing the boxes of the elements of 
A; (Hi) boxing the elements of B; then r = v. 

Take for instance 




and choosing for R the first row, we have A 
we obtain 



{a, b,c,d}, B = {q,/3,7, 5, e} and 




a bed 

H / 9 
T| j k I 



E 

h i 



m 



Note that here r represents a sum of 12 products of minors, whereas v stands for 
a sum of 30 products. 



4. Wronskian type determinants of orthogonal polynomials 

We are now ready to establish and prove the main result of this paper. 

Consider a linear operator T acting in the linear space of polynomials P and 
satisfying that deg(T(p)) = deg(p) — 1, for all polynomial p. We associate two 
sequences of polynomials (r n )„ and (s n ) n to this operator T. The sequence (r„)„ 
satisfies that 



(4.1) 



ro = 1 and the degree of r n is n, n > 0, 
T(r„) = r„_i, n > 0. 



This sequence of polynomials always exists and it is unique if we fix the values 
of r n , n > 1, at a given number xq. Indeed, fix a sequence of numbers (£ n )n>i 
and assume that we have constructed the polynomials Tj, j = 0, • • • , n, satisfying 
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rj(xo) — £j and (|4.1|) . Since T((x — xo)™ +1 ) is a polynomial of degree n, we can 
write T((x — xo) n+1 ) — X^=o a i r h wl th a n ^ 0. We then define 



r n +i[x) 



(x - x ) n+1 ~ £™ 1 otj(rj+i(x) - Cj+i) 



+ 1- 



From where we easily get that r„ +1 has degree n + 1, T(r„ +1 ) = r n and r n+1 (a;o) = 
The sequence (s„)„ are now defined recursively by so = 1 and 

n 

(4.2) Y,s j (x)r n - j (x) = 0. 

It is easy to see that if we write ^ r (x,t), fy s (x,t) for the (formal) generating 
functions of the sequences (r n ) n , (s„) n , respectively, 



(4.3) 



* r (x,i) = J2 r n(x)t n , *s(x,t) = ^S n (x)t n , 



n=0 



n=0 



we have \& r (a;, t) i S s (x, t) = 1. 

Theorem 4.1. Lei T be a linear operator acting in the linear space of polynomials 
satisfying deg(T(p)) = deg(p) — 1, for all polynomial p. We associate it the two 
sequences of polynomials (r„)„ and (s n ) n defined as above (see $4-l\ ) md |^.l?| ) ). 
Let (, ) be a bilinear form acting in the linear space of polynomials. For a given 
sequence of polynomials (ipi)i, tpi of degree i, we write fij = (rj,ipi), < 
With the operator T and the bilinear form (, } we associate the polynomials p n , q l n , 
n,i > 0. The polynomials p n , n > 0, are the pseudo orthogonal polynomials with 
respect to (, ) defined by Pq — 1 and for n > 1 



(4.4) 



Pn(x) = 



H'n-l H'n-l 



/4! 



,,n-l r 
rn ' n 

r n-l 



^r 1 



The polynomials q l n , n,i>0, are defined by 
(4.5) 

We then have 



3=0 



(4.6) det {r-\ Pm+J - 1 )) =1 = C n>m det U^n-i) 

where C n . m — (— 1)( 2 ) fl?=o ®"m+j is independent of x (and Q n is defined by (2.$fy ). 



Proof. The proof is similar to that of Leclerc for Theorem 11.11 in the Introduction 
(see [HI). 
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Let us introduce the (to + n) x (to + 'in — 2) matrix 



/ 1 ■■■ /C + „_! fm+„-i ■■■ tC+l-l r m+n -i T(r m+n -i) ••• X 1 " _1 (r- m +„_i) ^ 
1 ■■■ /C + „_ 2 |U m+ „_ 2 ■■■ Mm+nl2 r m+ „_ 2 T(r m+ „_ 2 ) ••• T"- 1 (r m+ ™-2) 



■■■ 1 flUl ^n« + l 

■■■ Mm 



T— X (r m ) 



\0 ■■■ Mo ■■■ M™ + "" 2 1 0-0 

The column vectors of this matrix will be denoted from left to right by: 

1,2,..- .n-l^V 1 -.- , Ai m+ "- 2 ,r,T(r),... .T"" 1 ^), 
Using the notation of Section 3, the left-hand side of (|4.6|) is written as 

.0 i ,7H— 1 v.1 Mo i'i O/(0 ,,m t 



,0 ,,m+n-2. 



[12 •••n- 1/i ---^ r] [12---n-2^ u ---^ m r] ■■■ [/r • • • /x""-— rj 

[12 ■••n- l A i°---^ m - 1 T(r)] [12 • • • n - 2^° • • • /x m T(r)] ■•■ [/x° • • • ^ m +"" 2 T(r) 



[1 • • • n - 1/ • • • fj, 



m — 1 rpn 1 



(r)] [1 • • • n - 2/i° • • • ^T"- 1 (r)] 



Now using the transformation of Theorem 13.11 with i? being the last row, this 
tableau is equal to 



2 
2 



n 
n 



n- 1 



M 1 A* 2 



T(r) 



A* 



,m+l 



..m+n — 3 ..m+n— 2 



m-2 / 



1 


2 


• n - 2 n- 1 


A* 




. . m — 1 . . rn 
fJ, ^ 


1 


2 


• n - 2 /i° 


A* 1 




m ,.rn+l 
(J, jJ, 


1 


> • 








^m-\-n — 3 ijTn+n— 2 


M° 


A* 1 • 


M »n-1 


r 


T(r) • 


•• T n ~ 2 (r) T"- 1 ^) 



Note that all boxes have been deleted in the second tableau because all the other 
permutations of the letters enclosed give rise to tableaux with two equal letters on 
some row, which are therefore equal to zero in view of the skew-symmetry of the 
determinant. Thus, we have rewritten the Wronskian type determinant of (I4.6P as 
a product of n determinants, where only the last one depends on the variable x. 
Explicitly, we have obtained that the left-hand side of (|4.6I) is equal to 



^m+n— 2^m+n — 3 ' ' ' ^m^(^) 
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where 



F'm+n-l 
11° 



m? 
m8 



Mrn+n-1 
Mm+n-2 



Ml 

Mo 



m-l 
H'm+n-l 

m-l 
fm+n-2 



Mi 

M™ 



' m+n— 1 
r m+n-2 



1 



T(r m+n -x) 

1 





T'^H^+n-l) 
rpn— 1 / 



m+n— 2 J 



and Q„ is defined by (|2.2p . Taking into account that T(r n ) = r n _i, we have 



e(x) 



Mm+n-2 



Ml 

m8 



Mm+n-1 
Mm+n-2 



Ml 
Mo 



m-l 
r-m+n- 

m-1 
H'm+n- 



Ml 

McT 



'"m+n— 1 
?"m+n-2 



ri 
1 



?"m+n-2 
?"m+n — 3 

1 





We now add to the first row of this determinant the linear combination of the next 
m + n — 1 rows multiplying the j-th row by Sj-\{x), j = 2, ■ • • ,m + n. (14.51) and 
(|4~2)) then give 



9(x) 



y?n+n~ 
f^m+n- 



11° 

Mi 



m8 



*?m+n— 1 



m! 
Mo 



^Vn+n— 1 

rn-1 
^m+n-2 



Ml 

McT 







n. 
1 



o 

'm+n— 3 

i 

o 



o 

r m -i 






In a similar way, we add to the second row of this determinant the linear combina- 
tion of the next m+n — 2 rows multiplying the j-th row by Sj-2(x), j = 3, • ■ • , m+n. 
(|43]l and (|4~2|) then give 



0(1) 



Im+n-l 
n° 

Hm+n-2 
IJ° 

P-m+n-3 



M? 

m!! 



^Vn+n— 1 
Qm+n-2 
Mm+n-3 



Ml 
Mo 



m-l 
"m+n— 1 

m-l 
ym+n-2 

m-l 
A'm+ra-S 



Ml 

M™ 






'm+n— 3 



r-i 
1 






'm+n— 4 

l 

o 



o 



Tin— 1 

o 





Proceeding in the same way m times we have 



Im+n- 



H-n- 



9m+n- 



Mn- 



/'{-! 



Mo 



And hence 



e(a;) = 



lm+n-1 



Mn-1 



'm+n- 1 "m+n-1 



^m+n— 1 



~m — 1 



This finishes the proof of the Theorem. 



□ 



12 



ANTONIO J. DURAN 



Theorem 11.21 in the Introduction is an easy corollary of Theorem 14. II just taking 
into account the expression (|2.3|) for the monic orthogonal polynomials p n with 
respect to a measure and that p n = £, n Pn- 



5. Specialization to some relevant operators 

For some relevant operators the estructure of the determinant in the right hand 
side of (|1.7|) is richer. Indeed, as we can see in Leclerc identity fj 1 .3|) . for the 
derivative the polynomials ipi, i > 0, can be taken so that this determinant has 
a Hankel structure. That also happens for the q-derivative. For the first order 
linear operator A there are also choices of the polynomials ipi, i > 0, so that this 
determinant has almost a Hankel structure. There also are other operators for which 
the polynomials tjji, i > 0, can be taken so that the determinant in the right hand 
side of (11.71) has a Toeplitz structure. The operators defined by T M (p„) = p n _i, 
where {p n ) n is a sequence of orthogonal polynomials with respect to any given 
measure /z, are among these operators. 

On the other hand, these relevant operators have associated distinguished se- 
quences of orthogonal polynomials for which the polynomials ipi, i > 0, can be 
taken so that the polynomials q^, n,i > 0, are again orthogonal. That is the case 
of the classical families of Laguerre and Jacobi for the derivative and the classical 
discrete families of Charlier, Meixner, Krawtchouk and Hahn for the first order 
difference operator. For these examples, the polynomials q l n belong, up to a change 
of variable, to the same family (but with different parameters). The case of the 
q-derivative is more involve: when one considers the q-classical families the poly- 
nomials q l n might not belong to the same family but to a relative one. For certain 
measures \x (for instance, Laguerre, ultraspherical, Charlier and Meixner measures) 
and the corresponding operator T M (defined as above), the polynomials ipi, i > 0, 
can also be taken so that the polynomials q l n , n,i > 0, belong, up to a change of 
variable, to the same family (but with different parameters). 

5.1. The derivative. In the case of the derivative, T — d/dx, we can take 

(5.1) r n (x) = 



ii 



(5.2) ^(x) = (x-x y, 

where xq is a fixed number. A simple computation then shows that 

(5.3) s n (x) = : , 

(see (|4.1[) and (|4.2I) for the definition of the polynomials r n , s n , n > 0.) If we write 
Mn — J{ x ~ Xo) n dn for the generalized moments of the measure /z, we then have 

(5.4) = J ipiTjdfj, = J {x -*f +3 dn = >±-L. 

From where it easily follows that q^ 1 — (x — xo)q l n — (n + l)q' l n+1 , and for x = 
q n = n\q^, where q n , q l n , n,i > 0, are the polynomials defined by (jl.2j) and (|1.6|) . 
respectively. This shows that the determinant in the right hand side of (jl.7| can be 
reduced to the Hankel determinant in the right hand side of Leclerc identity (|1.3p . 
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We now check that for the classical families of Laguerre and Jacobi, the number 
xq and the sequence of polynomials (ip n )n can be chosen such that the polynomi- 
als q l n belong, up to a change of variable, to the same family (but with different 
parameters) . 



5.1.1. Jacobi polynomials. For a, f3 S K, a + (3 ^ —1, —2, • • • , we use the standard 
definition of the Jacobi polynomials (P" )n m terms of the hypergeometric function 



(,5> ^w-(T)£ '"^::tr% >-'> 



j 



(that and the next formulas can be found in [6 , pp. 169-173). 

For a, j3,a + /3 ^ — 1, — 2, • • • , they are orthogonal with respect to a measure 
Ma,/3 = Ha,0{x)dx, which we normalize by taking 

^(x)ax-Z T (a + P + 2) ■ 

We write c Qji g for the constant in the right hand side of this identity. 
The measure ^t Q ,/3 is positive only when > —1, and then 

(5.6) fi a ,f;(x) = (1- x^il + xf, -l<x<l. 

We now show that by choosing xq — 1, the polynomials q l n , n,i > 0, are, up to 
a change of variable, again Jacobi polynomials (but with different parameters). 
Indeed, we have that 

o [, , v, (-2) n (a + l)„ 

J (a + p + 2) n 

Taking into account the definition of the polynomials q % n (| 1 . 6|) and (15 ,4[) , we have 
A(-2)^( a + l) l+J (l-x)"^ 



<1 



. = - (a + p + 2) i+j j\ (n-j)\ 



(-2)'(q + l) t v (-2V(q + i + l) 3 n\ 

The definition (|5.5|) for the Jacobi polynomials gives 

'n + a + 13 + i+l 



(1 - a .)»pa+/»+i+l,-»-/»-l f = ( ? 



7? 



(5.8) x A(-2p-(. + , + lW^ 

Comparing (|5.T[) and (|5.8p . we deduce that 

(a + /3 + 2) rl+4 V^^ 1 / 
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5.1.2. Laguerre polynomials. For a 6 R, we use the standard definition of the 
Laguerre polynomials (L") n 



(5.9) 



3=0 



n + a\ (—x) 3 
n-jj j\ 



(that and the next formulas can be found in [6], pp. 188-192). 

For a =/= — 1, — 2, they are orthogonal with respect to a measure \x a 
fj, a (x)dx, which we normalize by taking 



fj, a (x)dx = T(a + 1). 

The measure fi a is positive only when a > — 1, and then 
(5.10) fj, a (x) = x a e- x , 0<x. 

One can see that by choosing xq — 0, the polynomials <f n , n,i > 0, are, up to 
a change of variable, again Laguerre polynomials (but with different parameters). 
Indeed, proceeding as before for the Jacobi polynomials, we deduce that 



= (-l)T(a + i + l)L; 



-ct—i—n—1 



(-X). 



5.2. The first order difference operator A. For the first order difference op- 
erator, T = A, we take 



(5.11) r n (x) = 

A simple computation then shows that 



(5.12) 

Using (|2.5[) . we get 



s n (x) 



—x 
n 



i-l 



J 

i - I 



= E 

1=0 

Using now ()2.6|) . we find 

In = E ( 

n+i min(n,g) 

= 5>2 E 



J+l 



EE 

j=0 (=0 



rf ^ = E 

2=0 

d^=E 



x — j ^ /x 
I 



dfi 



x 
J + l 



1=0 



J 

i - 1 



J + l 



J 

i - z 



t + j-g; \g-jj 



—x 
n-j 



E 



9-' 



— ,T + Z 

n + i — g 



(5.13) 



E^ 

3=0 



o ( g + i\ j — x + 1 



n-g 
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This formula for <f n allows us to show that the determinant in the right hand side 
of (|1.7[) has almost a Hankel structure. Indeed, a simple computation using (|5.13|) 
gives that 

f I + m «n(«) = T^n+lfc - X )- 

« + 1 i + 1 

If we write q n — nlq^, it is easy to see that the determinant in the right hand side 
of (|1.7p reduces to 

■ det (q n+i+ j- 2 (x - j + 1)) ■" 



n! m (l + i) m ( m_1 )/ 2 ^'wtttj-A" j 1 v/i.jsi- 

For the classical discrete families of Charlier, Meixner, Krawtchouk and Hahn, 
the polynomials q^, n, i > 0, are actually polynomials of the same family (but 
with different parameters) in the variable —x + i. This will allow us to prove our 
conjectures in |4] on Casorati determinants for Charlier, Meixner and Krawtchouk 
orthogonal polynomials. 

5.2.1. Charlier polynomials. The monic Charlier polynomials are defined by 
(5.14) £ n( x )=p^- a ) n - j (^)(*)^ a ^°> 

(that and the next formulas can be found in [3], pp. 170-1; see also [13], pp, 247-9 
or [T5], ch. 2). The Charlier polynomials (c") n are orthogonal with respect to the 
discrete measure 



(5-15) Pa = £ ^S. 



The norm of the monic Charlier polynomials is given by 
(5.16) |K|| 2 = e<W\ 

We now show that the polynomials q l n , n,i > 0, are, up to a linear change of 
variable, again Charlier polynomials (but with parameter —a). Indeed, a simple 
computation gives //!■ = J n)dp a = e a a? Using (|5.13j) we get 



g+i 



g + i \ I —x + i 



:Q (g + i)\\ i J\n-g 

i\n\ ^ \n-gj \ n-g J 

(5.17) = e ^L d -*(-x + l ). 

i\n\ 

The identity (|1.4p in the Introduction of this paper, or equivalently, conjecture 1.1 
in [I] , now follows as a consequence of the identity (|1.7[) for the first order difference 
operator A, (|5.17[) and the duality 

(-1)" f-l) m 

L_J_^ m = L_J_^ n , n,m>0. 

a n a m 

Corollary 5.1. Conjecture 1.1 in [3] (see holds for n,k,m > and a £ 

R\{0}. * 
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5.2.2. Meixner polynomials. The monic Meixner polynomials are defined by 

<«*> w-^t^m* 

(we have taken a slightly different normalization from the one used in [3J, pp. 175-7, 
from where the next formulas can be easily derived; see also [T3J, pp, 234-7 or [IS] , 
eh. 2). For a ^ 0, 1 and c ^ 0, — 1, —2, • ■ ■ , they are always orthogonal with respect 
to a moment functional p a .c, which we normalize by taking (/0 OjC , 1) = T(c). For 
< |a| < 1 and c ^ 0, —1, —2, • • • , we have 

(5.19) Pg|C= (i- )cg g!££l±^ 8 . 

x=0 

The moment functional /9 aiC can be represented by a positive measure only when 
< a < 1 and c > 0. The norm of the monic Meixner polynomials is given by 

(5.20) IK'T= n!r(c+n). 
A simple computation gives 

o qJT( C + j) 
' (1 - ")•*'./• ' 

Proceeding as before for the Charlier polynomials, we can prove that the polyno- 
mials q % n , n,i > 0, are, up to a linear change of variable, again Meixner polynomials 
(but with different parameters): 

(5.21) q i(x) = ^ + i \ rh a n '- c - n - i+1 (-x + i). 

«!n!(l — ay 

Conjecture (3.7) of [4] now follows as a consequence of the identity (II. 7[) for the 
first order difference operator A, (|5.21j) and the duality 

(a — ^) n (a — 1 ) m 

a»(l+c) n _i " v ; a m (l + c) m _! mV h ~ 

Corollary 5.2. Conjecture (3.7) in [3] ZioWs /or n, k,m > and a ^ 0,1. 

5.3. Krawtchouk polynomials. For a ^ 0, —1, we write (fc^ ,Ar ) n for the sequence 
of monic Krawtchouk polynomials defined by 

(5 , 2 ) <»w-^g(-iy (i+ - ),( ;-" ) - , (;)C' 

For a ^ 0, — 1 and iV ^ 1,2, they are always orthogonal with respect to a 
(signed) measure p a .N which we normalize by taking {p a ,N, 1} = 1- When N is a 
positive integer and a > 0, the first N polynomials are orthogonal with respect to 
the positive Krawtchouk measure 

T(N) ^ a x 

(5 - 23) Pa > N = (i + a )»-i ^ r(J\T 



A simple computation gives 



(1 + a)* j! ' 
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Proceeding as before for the Charlier polynomials, we can prove that the polynomi- 
als q^, n, i > 0, are, up to a linear change of variable, again Krawtchouk polynomials 
(but with different parameters): 

Conjecture (5.3) of [4] now follows as a consequence of the identity (|1.7|) for the 
first order difference operator A, (j5.24j) and the duality 

(" + !)" h a,N ( v (a+l) m ia>N , s . n 

— — —7 k' (m) = 7- — k' (n), n,m>0. 

o n (2-JV) n _i " v ; a m (2-7V) m _! m w ' - 

Corollary 5.3. Conjecture (5.3) in [I] holds for n, k,m > and a / 0, — 1. 

5.3.1. Hahn polynomials. For a + c — N ^ —1, —2, • • • we write (h"' c ' N ) n for the 
sequence of monic Hahn polynomials defined by 

(5.25) /C c '» = E (l-^ + jWc + jOn-,- / »\ / *\ 

^ (n + a + c-N + j)^ 

(we have taken a slightly different normalization from the one used in |18j . pp. 
30-53, from where the next formulas can be easily derived; for more details see |13j . 
pp, 204-211). 

Assume that a + c — N + l,a — N + l,a + c, c^O, —1, —2, • • • . If, in addition, 
N is not a positive integer, then the Hahn polynomials are always orthogonal with 
respect to a moment functional p a .c.N, which we normalize by taking 

T(a + 1- N)T(c)r(a + c) 

\Pa,c,N; 1) — 7T7 1 7~\ \T\ ' 

1 (a + c + 1 — N) 

When N is a positive integer, the first N Hahn polynomials are orthogonal with 
respect to the Hahn measure 

(5-26) Pa , c , N = T(N) ^ r(N _ x)xl ^ 

(which it is positive when a > N — 1 and c > 0). The norm of the monic Hahn 
polynomials is given by 

(5 27) n^cJVip _ nl(N-n) n T(a- N + n + l)T(c: + n)T(c: + a + n) 

[ ' ' 11 " 11 (a- N + c + 2n)T(a- N + c + n)(a - N + c + n)*' 

A computation gives 

o = (N- jjjTfr - JV + l)r(c + j)r(c + a) 
^ T(c + a-N + j + l)j\ 

Proceeding as before for the Charlier polynomials, we can prove that the polyno- 
mials q l n , n,i > 0, are, up to a linear change of variable, again Hahn polynomials 
(but with different parameters): 

1n\ x ) — c a,c,N,n,in n { — X -f- IJ, 

where c atC ,N,n,i is a constant independent of x. 

This can be used to establish some symmetries for Casorati determinants of Hahn 
polynomials different to the ones pose in the Conjectures 5.1 in [3J. 
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5.4. The q-derivative. For a number q ^ 1, consider the q-derivative T = D q 
denned by 

( , 28 , w-*^- 

We can take 

(1 - q) n x n 



(5.29) r n (x) - . 

(q;q)n 

(5.30) ip i (x)=x i . 
A simple computation then shows that 

qfi)(l - q) n {-x) r 



(5.31) s„(ar) 



(q;q) n 

If we write /i„ = J x n dfi for the moments of the measure /i, we then have 



(5.32) = / ijjirjd^L 



(1 - qp'/i 



(q;q)j 



From where it easily follows that 



1 — q™ +1 

«» - q"^ = 1 _ <?n+V 

If we write q n — o° , we can see that the determinant in the right hand side 

of (|1.7[) is equal (up to a multiplicative constant independent of x) to the Hankel 
determinant det(q n+l+j ^ 2 )™ j= i ■ 

5.4.1. q-Laguerre polynomials. For a =/= —1, —2, • • • , we use the standard definition 
of the q-Laguerre polynomials (L"- q ) n 

(5.33) L« W = £ (q-^)* 

(q;q)n ^ (q + ;q)j(q;q)j 

(that and the next formulas can be found in [13], pp. 522-525). 

Using that the moments of the q-Laguerre polynomials are given by 

r(-a)r(a + l) (l-q)"+ 1 (q a+1 ;q)n . . 
Lq(-a) q na +{ 2 ) 

one can see that the polynomials q % n , n, i > 0, are, up to a change of variable, little 
q-Jacobi polynomials ([13j, p. 482): 

q l n {x) = c nti , a p n (-q n - 1 x;c\- n - a - i ,0\q), 

where c n ,i, a is a constant independent of x. 
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5.5. The operator T^{p n ) = p n -i- Let (p n )n be a sequence of orthogonal poly- 
nomials with respect to a positive measure /i, and consider the linear operator T n 
defined by linearity from T M (p n ) = p n -i- Notice that the determinant in the left 
hand side of (|1.7p has now a Toeplitz structure. 
We take 

r„ (x) = p n (x) , ipi (x) = pi (x) . 

Define finally the polynomials s n , n > 0, by (I4.2[) . 
We then get 

\\Pi\\ 2 , if J = i, 



a*: 



= / ^iTjdn = \ PiPjdfi = 



if j ^ i 



This gives 
(5.34) 



3=0 



-j = lift 



Hence the determinant in the right hand side of (jl.7p has the Toeplitz structure 



s n (x) s n -i(x) s n - 2 (x) 

s n+ i(x) s n (x) s n -x(x) 

^n+m- 1(^) ^n+m- 2 (^) ^n+m— 3 (^) 



^n— m+2 (*^) ^n— m+1 j 
$n— m+3 ^n— m+2 

s„+i(ic) S„(l) 



For some families of orthogonal polynomials, the polynomials q l n (|5.34[) . n, i > 0, 
are, up to change of variable, actually polynomials of the same family (but with 
different parameters). Among these families are Laguerre, ultraspherical, Charlier 
and Meixner polynomials. 

5.5.1. Charlier polynomials. Since A(c£(x)) = nc°_ 1 (x), the determinants in (|1.7p 
for the operator T M are essentially the same that for the first order operator A (see 
Section |5~2~T|) . 

5.5.2. Laguerre polynomials. The generating function for the Laguerre polynomials 
dSH) is 

OO 

Using (|4.3|) we get that s n (x) = L~ a ~ 2 (—x), and, hence, taking into account (|5.34p 
we deduce that also q l n (x) = \\Lf \\ 2 L~™7 2 {-x). 

5.5.3. Ultraspherical polynomials. Consider the ultraspherical (o Gegenbauer) poly- 
nomials (C^)n, A ^ 0, defined by the generating function 

oo 

(5.35) (l-2xz + z 2 )- x = J2c^(x)z n . 

n=0 

If A = —k, with k a positive integer, we have C~ k — 0, n > 2k. It turns out that the 
finite sequence of polynomials C~ k , n = 0, • ■ ■ , k, satisfies a three term recurrence 
relation and hence they are orthogonal with respect to a (signed) measure. 
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For A = 0, the ultraspherical polynomials reduce to the Chebyshev polynomials 
of the first kind defined by the generating function 



1 — xz \ -\ 



1 — 2xz + z 2 ' 

n— 



The ultraspherical family (C„) n is the particular case of Jacobi polynomials for 
a = /3 = \- 1/2. 

Using (|4.3p we get that s„(a;) = C~ A (:r), A ^ and, hence, taking into account 
(|5.34[) we deduce that also q % n = WC^^C^^x). In particular, for A = k with k a 
positive integer, we get s n = 0, n > 2k + 1. 

For the Chebyshev polynomials of the first kind, it is easy to see that so = 1, 
s\(x) = —x and s n (x) — x n ~ 2 (l — x 2 ), n > 2. And for the Chebyshev polynomials 
of the second kind, we have sq = 1, Si(x) = —2x, S2{x) = 1 and s n {x) = 0, n > 3. 

As a consequence det ) = for n > 3 and m > 2. Theorem 11.21 then 

implies that for (p n ) n the Chebyshev polynomials of the first or second kind, also 
det [p m +3 -i( x )))" j—i = f° r n > 3 and m > 2 (where we take p n — for n < 0). 
This is, anyway, an easy consequence of the fact that Chebyshev polynomials of the 
first and second kind have, essentially, constant recurrence coefficients. 

5.5.4. Meixner polynomials. For a ^ 0, 1, consider the Meixner polynomials defined 
by the generating function 



n=0 

A straightforward computation shows that 

(a — l) n 

a.c/ \ V / ; a,c/„\ 

m n W = ^~r m n (X), 

a n n\ 

where mj c (i) is the n-th monic Meixner polynomial defined by (|5.18p . 

Using (|4.3p we get that s n (x) = ro£'~ c (— x), and, hence, taking into account 
(15.341) we deduce that also q l n — ||m™' c || 2 m™'~j C (— a;)- 

6. SELBERG TYPE INTEGRALS AND SUMS 

Write x for the multivariable x = (x\, ■ ■ ■ ,x m ) and A(x) for the Vandermonde 
determinant: 

(6.1) A(x)= H ( Xi - Xj ). 

l<i<j<m 

Given a measure consider the multiple integral 

(6.2) |...|s(x)A 2 ^(x)^(x), 

where d/i(x) — dfj,(x\) ■ ■ ■ d/i(x m ) and s(x) is a symmetric polynomial in the variable 

1 ' ' ' ? Xjji . 

When (i is the Jacobi measure dfi — x a_1 (l — x) l3 ~ 1 dx, a, ft > 0, and s(x) = 1, 
the integral f|6.2|) is the celebrated Selberg integral ([IS], see also [7]): 

'■' '■' "' '"^r(a + i7)r(/3 + i7)r(i + 0- + i)7) 



A^( X )Y[x^-\i-x j f- i d X = ii 



u r(a + /3 + (m+i-l) 7 )r(l+ 7 ) 
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As we wrote in the Introduction, other choices of s(x) give well-known extensions 
of the Selberg integral. For instance, s(x) = IIJLiC^j — u ) i s Aomoto's integral 
PP, and s(x) = J\{x\, ■ ■ ■ , x m , I/7), where 7 € N and Ja is a Jack polynomial, is 
Kadell's integral [TO]. For s(x) = 1 and fx the Hermite measure dfi = e~ x / 2 dx, the 
integral (|6.2p is the also celebrated Mehta integral. 

In this Section, we prove Theorem 11.31 in the Introduction. That is, we show 
that for certain operators T the determinant in the right hand side of (|1.7[) can be 
rewritten as an integral of the form (|6.2p for 7 = 1 and s(x) = njli r n{ x j)- Hence, 
the determinant in the left hand side of (|1.7p gives the value of this integral. 

This will allow us to obtain the value of some new Selberg type sums. 

The main idea is to use a variant of the Heine's integral formula ([5]) for a Hankel 
matrix: If fj, is a measure on the real line and we write fj,j — J x J dfi, j > 0, for its 
moments, then 



det((Mi +i - 2 )^ = i) = J ■■■J A 2 (x)d/i(x). 
We need the following lemma 

Lemma 6.1. Let fi be a measure in the real line and consider a sequence of poly- 
nomials (ip n ) n with ip n of degree n and leading coefficient equal to v n . Then, for 
any symmetric polynomial s(x) we have 
(6.3) 

r r m n^T 1 v- r r 

j ■■■j s(x)A(x)J[^ 1 (x j )dfi(x) = 3 j ■j S (x)A 2 (x)d M (x). 

Proof. For each permutation a of {l,--- , m}, we perform in the integral in the 
left hand side of (|6.3[) the change of variable = x a u\ . Since s is a symmetric 
polynomial and A(y) = (— l) e ^ <T ^A(x), where e(cr) denotes the signature of the 
permutation er, we get 

/P TTX p p TTL 

■■■ S (x)A(x)n^-ife)^(x) = (-l) eM /••• / s(x)A(x)J]^_ 1 (o ;ct(j) )^(x). 
J j= i J J j= i 

Summing up over all permutations we find 

J ■•• | S (x)A(x)n^-ife)^(x) = ^ / ••• / S (x)A(x)^(-l) e Wn^-i(^ W )^(x) 

S (x)A(x)det((^-_ 1 ( a ; 4 ))™ =1 )d/i(x). 

The latter determinant can be easily transformed into a Vandermonde determinant 
by elementary column operations (see [J2], Prop. 1) to get det((V>j-i (2^))™ =1 ) — 

(-l)C?) (njl^o 1 Vj) A ( x )- The P r001 is now finished. □ 



11 ' 



Given a number u and an operator T, we associate to T the polynomials r„, s 
n > 0, as in (|4.1I) . (|4.2I) . In this case we take r n (u) = 0, n > 1 (these conditions 
define r n , n > 0, uniquely). We write r n — r rhU to stress the dependence of r n on 
u. As an easy consequence we have that also s n {u) = 0, n > 1, (see (|4.2I) ). This 
gives that q r n (u) = tf n = J ipir n d^L, where q l n are the polynomials given by (|1.6p . 
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Hence the determinant in the right hand side of (|1.7|) evaluated at x = u can be 
rewritten as 



det («i-K-iM). j=1 = det yj ^j-ii^rn+i-iA^Wixj) 

(6.4) = J ... J ( JJ^-_ 1 (a; i )J det (r n+ i_i, u (o; i ))^. =1 d/x(x). 

We are now ready to prove Theorem 11.31 in the Introduction. 
Proof. We first prove that 

(6.5) det U- + U{u)) m = [ ■■ / ft ' W rn,»fe)W, 

3=1 

where u n , cr„ are the leading coefficients of the polynomials ip n and r n , respectively. 
Since r n<u (x) = a n Hj =1 (x - %•), we get 

m 

(6.6) det (r n+J _i(x J ))™ =1 = det (^-lO^-))™ =1 r ra , M (a;j) 

3=1 

where /ij(ar) = nj=ri+i( a; ~ a i)- The determinant in the right hand side of ()6.6[) 
can be easily transformed into a Vandermonde determinant by elementary column 

operations (see [12], Prop. 1) to get det (/ij_i(a;j))™. =1 = (— l)( 2 )A(x). Hence 
(jO|) and |JO]) give 



... / H^fe) A(x)U^-i(*i)dMx). 

J \ i=i °" / j=i 

(|6.5[) now follows easily using Lemma [B~T1 The proof of the Theorem can be finished 
by using Theorem 1 1 . 2 1 and the formula (|2.4[) for Q„ in terms of the norm of the monic 
orthogonal polynomials. □ 

The three relevant operators T — d/dx, T — A and T = D q satisfies the hypoth- 
esis of Theorem 11.31 with 

_ (x - u) n _ fx - u\ (1 - q) n x n (u/x;q) n 

r n,u(x) — : , i nu (X) — i i, ana r nu (x) — , 

n\ \ n J (q;q)n 

respectively. Hence, we have 

Corollary 6.2. Let fi be a measure and write {p n ) n for the sequence of monic 
orthogonal polynomials with respect to fi. Then 

/ • • • / A 2 (x) H( Xj - u)"d/x(x) = C d , hm det (p^_ 1 (u)) i =i , 
3=1 2J 

[ ■ J A 2 (x)f[(xj -u-n + l) n dfi(x) =C A ,n, m det (A i_1 (p m+j _i(u)))" j=1 , 

' ' 3=1 

J ■J A 2 (x)]^[x™(u/a; 3 ;q)„d At (x) =Cr.,,„, m det ( J D;- 1 (p m+J _ 1 (w)))™ j=1 , 
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where 

n;=oj! 

(-ir»m!(l-q)(S) nr^ife-f 
n j= o w; q)j 

respective/?/. 

When /i is one of the the classical weights of Jacobi, Laguerre or Hermite, and 
u = 0, the determinant in the right hand side of the first identity of the previous 
Corollary can be easily computed using standard determinant techniques. In par- 
ticular, for \x the Jacobi weight and u = 0, we recover the Selberg integral for 7 = 1. 
For fi the Jacobi weight and n — 1 we recover Aomoto's result for (|6.2[) when 7=1 
(PQ). In a similar way, for the classical discrete polynomials we get a deduction for 
Askey conjectures 5, 6 and 7 in [5] when 7 = 1. Similar results can be obtained for 
some of the q-classical weights and Askey conjectures for q-Selberg integrals. 

Theorem 11.31 can be rewritten only in terms of the polynomials r n . 

Corollary 6.3. Let (a n ) n , (o~ n ) n be two sequences of numbers with Cq = 1. Write 
r n , n > 0, for the polynomials of degree n defined by ro = 1, r n (x) = o~ n Ylj=i( x ~ 
a j), n > 1- Consider a measure fi and its associated sequence {p n ) n of monic 
orthogonal polynomials. For each n > 0, we can write p n (x) — X)J=o r j( x )- 
Then 



(6.7) J ■ J A 2 (x) l[ r n ( Xj )dn(x) = C T ^ m det (Cf 1 



where 



C Ttn , m = (-ir n m\o-™ J] \\ Pj \\ 2 llo- 

3=0 3=0 



n-1 



Proof. Define the operator T by linearity from T(r n ) = r n -i, n > 1, and T(l) — 
0. Then T l (p n )(x) = J2™ =0 0™rj-i(x), and since r n {a\) = 0, n > 1, we get 
= ^ is enough now to apply Theorem 11.31 □ 

Here it is an illustrative example of how to use this Corollary to get a Selberg 
type sum for Racah polynomials. 

6.1. Racah polynomials. Let N,(3,j,5 be real numbers satisfying that N is a 
positive integer and 5, 7, — (3 + 7 7^ — 1, • • • , —N, [3 + 8 ^ — 1, • • • , — 2N — 1. Define 
the numbers o~ n and the polynomials r„,r„, < n < N, by 

(-1)« 



(-JV) n (/3 + 5 + l) n (7 + l)„n!' 

n-l 

r n (a:) = (-1)™ J] (as - ifo + d + 1 + i)), 

i=0 



r„(x) = ct„(-1)"t„(:e). 
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Write X(x) — x(x + 7 + 5 + 1) and for < n < N consider the monic Racah 

polynomials (see [13], pp. 190-196) denned by 

(6.8) 



They are orthogonal with respect to the discrete measure 

>p (-N) x (f3 + 6 + l) a (7 + 1)^(7 + S + l) a (( 7 + 5 + 3)/2) a 

^ o N + 7 + S + 2) x (- j 8 + 7 + l) x (( 7 + 5 + l)/2) x (5 + l) B ar! A(a;) ' 

The norm of the n-th monic Racah polynomial is 

l]N N,^M2 = (-l)"M(/3 - 7 - JV) n (-<S - AQ»(/3 + l) n 
11 n 11 «T„C8-JV+l)2„C8-JV + n) n 

where AZ - { _ (3+7+ i )jv((5+ i )jv • 

A simple computation shows that r„(A(x)) = (— x) n (x+ , y + 6 + l) n . Hence, (|6.8[) 
can then be rewritten as 

With the notation of Corollary 16. 31 we have 

fl n_ (-^"(-^-(n + ^-JV),- 
3 n!cr„(/3 + n-7V)„ 

Using standard determinant techniques we get 

i! 



dot 



(>:":' ')" =(-i) mn n 



J= o ( m + j) !<7 m+i + w + j - iV) m 
Hence Corollary 16.31 gives the identity 

(-N,p + 6 + l,j+l,j + 6 + 3)/2) Wj (j + 5 + l) Xj+n 



E A2 (^)I1 



^ (N + 7 + S + 2, -0 + 7 + 1, (7 + S + l)/2, 5 + l) x . ( Xj - n)\ 

= , l)m n M m TT (J + W -l-N-5-N,p + 1U-N, (3 + S+l,J + l) n+j 
[ ' f} (f3-N + l) 23 (/3-N + j) 3 (l3-N + m + j) n 

where x A = (A(xi), • • • ,X(x m )), and (ui, • • • ,u k )j = (ui)j ■ ■ ■ (u k )j. 

Similar versions can be get for the other Racah weights (see [15] . pp. 190). 



7. Constant term identities 

In this Section 7, we show that for certain operators T and certain families of 
orthogonal polynomials (p n ) n , one (or both) of the determinants in (|1.7[) can be 
rewritten as the constant term of certain multivariate Laurent expansions. Some 
of them are generalizations of Morris identity (|1.11[) for k = 1. We next consider 
illustrative examples of this for Meixner, Charlier and ultraspherical polynomials. 
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7.1. Meixner polynomials. For a ^ 0, 1, consider as in Section l5.5.4l the Meixner 
polynomials (m°' c ) n defined by the generating function 

oo 

(7.1) (l - ~Y (1 - z)- x - c = m% c {x)z n . 

a n=0 

Summing up rows we have 

(7.2) det (A^K^-r (*)))" j=1 = det (m^O* + * - 1))^ =1 . 

Using Theorem 11.21 and taking into account the computations in Section 15.2.21 we 

get 

(7.3) 

det (a-v^ViW)) ■ J= i = c "-™ det (< + 7-T"^' +3 (-^+i - 1))". =1 - 

where 

_ (-i)" m+ (?)+(S)(i- )(5) 

^ n , m TnA 7mA * 

The generating function (|7.1j) gives for the Meixner polynomials the integral repre- 
sentation 

i r (i _ z) x (i _ z'l-a-c 
(7.4) <' C W = ^J C [ - Q M +1 ^ 

where the function /^(z) = (l — (1 — 2;) _a:_c is taken so that it is analytic at 
z = and / x (0) = 1, and C is a circle with center at z = 0, positively oriented 
and such that f x (z) is analytic in a neighborhood of C and in its interior. Using 
this integral representation, we can rewrite the determinant in the left hand side of 
(|7.3[) in terms of the constant term in the Laurent expansion of certain multivariable 
analytic function. 

Indeed, using (17. 4j) the determinant in the right hand side of (|7.2|) can be rewritten 

as 

■i-r.-ftofe-*)"'"^'-'-"')" *, 

V z j 3 /,,, i 



(2ni) n 

where dz = dz\ ■ ■ ■ dz n . This multiple integral can be computed as follows 



/ i _ £i — (i _ z y 
' V z 3 J,.,, : 



, i-1 ' 



= (-i)(s) / ... / fr . TT (1 ~ a)(Zt ~" j) , 



f [fr (l^jfT TT {z ._ z . ?dz 



where we have used the version of Lemma 16.11 for a contour integral 
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The Residue Theorem then gives that the last integral is equal to 

n /j _ zj_Y 

(27ri)™C.T.| z -oTT " — r — r TT 0* - z,) 2 , 

j — 1 \ 3' j l<«<J<n 

where if F(z) is a multivariable meromorphic function F(z) at = • • • = z n = 0, 
we denote by C.T.| z= oi 7 '(z) the constant term in the Laurent expansion at z\ = 
■■■ = Z n =0 of F(z). 

Using (|7.3|) . we finally get 

C.T.L-oTT " — 3— r TT (Zi-Zif 



LJ./-I _ _ \x+c+n-l m+n-L 
j = l^ Z 3> Z \<i<j<n 

(_l)(?)+(S)+»"» n ! a (?) det {rn^-^-'-^i-x + j - 1)) 



7.2. Charlier polynomials. Consider the Charlier polynomials (c°)„ defined by 
the generating function 

oo 

(7.5) {l + zfe-» = J2<(*)z n - 

n=0 

A straightforward computation shows that 

<(*) = c^(x)/n!, 

where c°(x) is the n-th monic Charlier polynomial defined by (|5.14[) . 

Proceedings as in the previous Section (using Theorem 11.21 and the computations 
in Section [5. 2 we get 

(7.6) det {c a m+3 ^{x + i - l));; j=1 = (-1)™ det (c^^-x + j 1))™ =1 . 

The generating function (17.5[) gives for the Charlier polynomials the integral repre- 
sentation 

(") ^-si^^*' 

where the function f x (z) — (1 + z) x e~ az is taken so that it is analytic at z = and 
/x(0) = 1, and C is a circle with center at z = 0, positively oriented and such that 
/x(^) is analytic in a neighborhood of C and its interior. 

Using this integral representation and proceeding as in the previous Section, we 
can rewrite the determinant in the left hand side of (|7.6p in terms of the constant 
term in the Laurent expansion of certain multivariable analytic function. In doing 
that, we get 

c.T.| z=0 n e -^^±gf n (^-^) 2 =(-i)® + " m -!dot(c-_ 1 (^ +J -i ) )^ 1 

j=l Z j l<i<j<n 

7.3. Ultraspherical polynomials I. Consider, as in Section 15.5. 3[ the ultras- 
pherical (o Gegenbauer) polynomials (C„)„, A ^ 0, — 1, — 2, • • • , defined by the 
generating function 



(7.8) (l-2xz + z 2 )- X = J2c^ 



,. , x > z 

n=0 
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Consider the operator T{C X ) = C x _ x . Using Theorem 11.21 and the computations 
in Section 15.51 and 15.5.31 we get 

(7.9) detCC^Oz))^! = (-ir-det^-^))^!- 

The generating function (|7.8|) gives for the ultraspherical polynomials (A ^ 0) the 
integral representation 

where the function f x (z) = (1 — xz + z 2 )~ x is taken so that it is analytic at z = 
and /a;(0) = 1, and C is a circle with center at z = 0, positively oriented and such 
that f x (z) is analytic in a neighborhood of C and in its interior. Notice that if we 
write C x (x) = for n < 0, the integral representation (|7.10l) also holds for n < 0. 

Using this integral representation and proceeding as in the previous Section, we 
can rewrite the determinant in the left hand side of (|7.9|) in terms of the constant 
term in the Laurent expansion of certain multivariable analytic function. In doing 
that, we get for A ^ 0, —1, —2, • • • 
(7.11) 

™ (i - 2xzj + z 2 y x „ „ fn\, > 

j = l Z 3 \<i<j<n 

Actually, this formula holds for A ^ 0. Indeed, for A = — k, with k a positive integer, 
the formula (17.8[) allows to define the polynomials C~ k , n > 0. For n > 2k we have 
C~ k = 0. It is not difficult to see that for fixed n, m, the determinant in the right 
hand side of (|7.11|) is a polynomial in A vanishing at A = 0. Since the left hand side 
of (|7.11[) is a continuous function in A and they are equal for A ^ 0, — 1, — 2, • • • , 
they are actually equal for all values of A ^ 0. For A = 0, the left hand side of 
(17.111) is then (this happens because the Chebyshev polynomials of the first kind 
do not satisfy the generating function f|7.8[) ^1 . 
For x — 0, the determinant in the right hand side of (|7. 1 1|) can be explicitly 
computed. In doing that, one has 

n (1 4- z 2 )~ x 

c.T.i^n 1 -^!- n (*-*) a 

j—l J l<i<j<n 

!0, if nm is odd, 

Lo(n, to, A), if nm is even and n < m, 

w(m, n, —A), if nm is even and rn < n, 

where for n < m we define 

n—l r ■ /ol I I^S^] -1 

W (n,m,A) = (-l)»m II (A+ J ) min( "' m - 2 -'' ) (A^-l) max( "- 2 '''- 2 ' 0) . 

3=0 y 2 \- j=0 

7.4. Ultraspherical polynomials II. If we consider the derivative instead of the 
operator T considered in the previous Section, the integral representation (|7.10[) 
leads to other interesting constant term identity. 
For ultraspherical polynomials we have the formula 



(C x )'{x) = 2\C x t\{x). 
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This gives 

(7.12) det ((^-i(*)) (i ~% =1 = 2 (S) (nW-O dct (^WIJh ■ 

Using Theorem ll.2l and the computations in Section l5.1.1l we get for A ^ —1/2, — 1 ; — 3/2, —2, 
(7.13) 

ft _ „-in-M-l p2A+j-l,-n-i-A+l/2 / 2+3 \ ' 
I 1 ^ r n+i-l ^i=T J 



((c4 +J --i(«)) (i - 1} ) B 



dot (JC; ■") =r„.,„d.-t 

where 



(2A + 2j - l)„ +l _ 3 



c n , m .(-ir(™) 2 (T)0|r 1(A) 



nr=7(2A), 

Using the integral representation (17.101) . and proceeding as in the previous Section 
we can write the determinant in the left hand side of (|7.13[) in terms of the constant 
term in the Laurent expansion of certain multivariable analytic function. In doing 
that, we get for A ^ -1/2, -1, -3/2, -2, • • •, 

™ (1 - 2xzj + z 2 A- x - n+1 „ o/l 

3=1 3 l<i<j<n V 1 3 



rn — 1 



(-l) nm +(?)n!2 nTO+ (») ( [| ] det 



f 1 _ ^"-H-l p2A+2-l,-n-i-A+l/2 / x+3 \ 
K X X ) r t»+i-l \x-lj 



For a; = 1, the determinant in the right hand side of the above identity can be 
explicitly computed to get 

n (1 v \-2A-2n+2 / 1 

c.T.i^oii 12 — ^ — n (*-*) a (i-— 

J - J - z' - L - L V ZiZj 

j = l 3 l<i<j<n v 1 J 

IJ !2 2 ™+(™) fnr=Ti 1 ( 2A + 2.? - 1; 



" 1 ^ j!(A)„ +i (A + i + l/2) r 



(n£ / i 2l (2A + n + 2 J -l) [ _ i]+[I2 + i] ) /i (n + j)!(2A) j (2A + 2i + l)„- 
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